Entanglement, Non-linear Dynamics, and the Heisenberg Limit 
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We show that the quantum Fisher information provides a sufficient condition to recognize multi- 
particle entanglement in a A/" qubit state. The same criterion gives a necessary and sufficient 
condition for sub shot-noise phase sensitivity in the estimation of a collective rotation angle 0. The 
analysis therefore singles out the class of entangled states which are useful to overcome classical 
phase sensitivity in metrology and sensors. We finally study the creation of useful entangled states 
by the non-Unear dynamical evolution of two decoupled Bose- Einstein condensates or trapped ions. 
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Introduction. The ability to create and manipulate en- 
tangled states of many-particle systems is a far-reaching 
possibility of quantum mechanics. Several efforts have 
been devoted, in the last few years, to exploit entangle- 
ment to design new technologies for secure communica- 
tion, metrology and fast computation or to unveil foun- 
dational problems of quantum mechanics. From the ex- 
perimental point of view, trapped Bose-Einstein conden- 
sates (BECs) [H, [2^, cold/thermal atoms [s^ and trapped 
ions |j| are important candidates for the creation of large 
scale quantum entanglement. It is important to empha- 
size, however, that not all entangled states are equally 
useful for developing protocols that outperform classi- 
cal operations. Generally speaking, current measures of 
entanglement mostly focus on the algebraic separability 
properties of quantum states. This notion should be ex- 
tended for quantum technological applications, where it 
is essential to classify entanglement on the basis of some 
additional physical/algebraic properties required by the 
specific task. These attributes are crucially related with 
non-separability, but are not necessarily possessed by all 
entangled states. 

In this Letter, we develop a general framework to study 
the interplay between entanglement and phase estimation 
in metrology and quantum sensors [5]. A quantum state 
Pinp must necessarily be entangled in order to be useful 
for estimating a phase shift with a sensitivity be- 
yond the shot-noise, which is the maximum limit attain- 
able with separable states. Nevertheless not all entangled 
states can perform better than separable states. Here we 
introduce a new criterion, on a generic pinp, which is suffi- 
cient to recognize multi-particle entanglement and is nec- 
essary and sufficient for sub shot-noise phase estimation 
sensitivity. We separate entangled states in two classes 
on the basis of an additional geometrical (or kinetic, see 
below) property in the Hilbert space. Our analysis uses 
basic tools of parameter estimation theory and provides a 
simple and experimentally measurable condition, Eq.([3]), 
which extends other criteria discussed in the literature 
based on the concept of spin squeezing [l]. We will show. 
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with an example experimentally achievable with dilute 
BECs and trapped ions, how non-linearity can generate 
a class of states which are entangled, useful for sub shot- 
noise interferometry, but not spin-squeezed. 

A state of N particles in two modes {N qubits) is sep 
arable (non-entangled) when it can be written as 



(1) 



where pk > 0, ^^.Pk = 1 and pl.^ is the density matrix 
for the ith particle. How to recognize entangled states? 
Let us introduce the "fictitious" angular momentum op- 
erator, J = '}2i=i ^^^^ where cr^^^ is a Pauli matrix op- 
erating on the /th particle. According to the current 
literature, if a state pinp satisfies the inequality 



iV(AJ„-3)2 



< 1, 



(2) 



then is particle-entangled 0, 0, 0] and spin squeezed [H, 
0, |l^ along the direction na, being ni, n2 and n^, three 
mutually orthogonal unit vectors and . = J - Hi. 

Here we introduce a different sufficient condition for 
particle-entanglement : 



N 



^Q[pinp, Jn] 



< 1, 



(3) 



where i^Q[pinp, J^l = MAR)'^ is the quantum Fisher in- 
formation (QFI) [lH, [ll, [ll, \U] and n is an arbitrary 
direction. The Hermitean operator R is the solution of 
the equation {R^ pinp} = ^[^n^ Anp] [15]. It is possible to 
demonstrate that ^ ^'^ [16]. Therefore, Eq.([3]) rec- 
ognizes a class of states which are entangled, < 1 
and not spin-squeezed, > 1 as, for instance, the max- 
imally entangled state ^]. Notice that, for a pure state. 



Pinp = |V^inp)(V^inp|, WC haVC Fqlpinp, Jfi] = ^^Jfi) [Ul 



and the sufficient condition for multiparticle entangle- 
ment, Eq.([3]), assumes the appealing form 



Xps 



N 



4(AJb)2 



< 1. 



(4) 



The QFI is naturally related to the problem of phase 
estimation. Generally speaking, an interferometer is 
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quantum mechanically described as a collective, linear, 
rotation of the input state by an angle 6: Pout{0) = 
e^^'^^Pinpg"*^'^^. The goal is to estimate 6 with a sen- 
sitivity overcoming the shot-noise limit AOsn = 
For instance, in Mach-Zehnder (MZ) interferometry, is 
a relative phase shift among the two arms of the inter- 
ferometer, and the rotation is about the n = y axis. 

For an arbitrary interferometer and phase estimation 
strategy, the phase sensitivity is limited by a fundamen- 
tal bound, the Quantum Cramer- Rao (QCR) [lH, which 
only depends on the specific choice of the input state, 

A^QCR = , ' „ = (5) 

A comparison with Eq.Q reveals that Eq.(|3]) is not only 
a sufficient condition for particle-entanglement, as al- 
ready discussed, but also a necessary and sufficient con- 
dition for sub shot-noise phase estimation. This is a main 
result of this work: x < 1 provides the class of entangled 
states which are useful for sub shot-noise sensitivity. In 
other words, chosen pinp and J • n, if the corresponding 
value of the QFI is such that x < 1, then the state is 
entangled and, if used as input of an interferometer re- 
alizing the unitary transformation e"*^'^^, it provides a 
phase estimation sensitivity higher than any interferome- 
ter using classical (separable) states. On the other hand, 
the class of entangled states for which x ^ 1 cannot pro- 
vide a sensitivity higher than the classical shot-noise. 

The QFI, which links Eqs.© and (O, has a simple 
interpretation as square of a "statistical speed", = 
^Q[Pinp,^n] = (d/(^)/d^)^. This corresponds to the rate 
of change of the absolute statistical distance l{6) among 
two pure states in the Hilbert space (or in the space of 
density operators for gene ral mixtures) along the path 
parametrized by [ll', 1^. The absolute statistical dis- 
tance is the maximum number of distinguishable states 
along the path parametrized by 6>, optimized over all pos- 
sible generalized quantum measurements. According to 
Eq.(|3]), useful entanglement corresponds to high speed, 
I'^fI > I'^crl, being \vcr\ = \fN a critical velocity that can- 
not be overcame by separable states Eq.([T]). The maxi- 
mum speed (strongest entanglement) is |^;max| = N and 
therefore the fundamental (Heisenberg) limit in phase 
sensitivity is A^hl = ^/N . Physically, this simply means 
that, under the action of some unitary evolution, use- 
ful entangled states evolve (become distinguishable) more 
rapidly than any separable state. 

Entanglement. Let us introduce the inequalities 

being Mk{e) = Tr[M^ 

Pout]? and M an arbitrary observ- 
able [l7]. The right-hand side of Eq.([6]) allows us to 
demonstrate Eq.(|3]) by showing that FQ[psep, Jn] ^ ^ for 
any arbitrary unit vector ft in the pseudo angular mo- 
mentum space. First, notice that, for separable states. 



Pk = Pk^ ^ Pk^ ^ ... ^ pi^\ we have 4(A4)2 = N - 
4Ez=i6l'V < Combining this result with Eq.® 
and the convexity of the QFI [16] (i.e. for an arbitrary 
mixture p = Y^j^pupki FqIp^Ju] < Y.kPkFq[pk, Jn]) we 
obtain that FQ[psep,^n] ^ where the equality sign 
can be saturated only with pure states. Moreover, since 
4(A4)2 < 4(J|) < N\ we obtain FQ[pinp, 4] < N\ 
Then, from Eq.Q, follows that A^hl is the highest pos- 
sible phase sensitivity. 

Using the left-hand side of Eq.® we now demonstrate 
that X ^ foi" arbitrary pinp. We consider, without 
loss of generality, a direction n = n2 such that (Jns) = 0. 
By choosing M = Jfi^ — (Jfts) in Eq.(|6]), we obtain that 
FQlpinpJn] > (dMi/d^)VM2 = TV/^^ Then, Eq.© 
shows that X ^ the class of states satisfying x < 1 is 
wider and includes the class of states deffned by Eq.([2]). 

Non-linear dynamics. We now discuss the connection 
between non-linear dynamics, entanglement and spin- 
squeezing. We consider a coherent spin state, \j^j)ni = 

EtL, ^V^(^|j,M>n3 E 0, with j = N/2. This 

state is separable (x^ = 1) and we quest the possibility 
to strongly entangle the particles by the non-linear evo- 
lution e ^3. A direct calculation of Eqs.([2]) and ([3|) 
with n = n2 (where the expectation values are computed 
over the state \ip{r)) = e '^^ ^^\jij)ni) gives 

e2 = (cosr)-2(^-i), (7) 

= 2I[{N +1)-{N- 1)(cos2t)^-2] . (8) 

Notice that, > 1, while < 1 for all values of 
r: the state |'0(t)) is not spin-squeezed but still (use- 
fully) entangled. A comparison between Eq.([7j) and 
Eq.® is presented in Fig.([Tla) for ^ 1. We em- 
phasize two time scales in the dynamical evolution of 
X^: for < r < I/VTV, X^ decreases from 1 to 
for l/>/]V < r < 7r/2 — I/a/TV, it reaches the plateau 
^2 _ 2/N. The dynamics are periodic with period 
T = 7r/2 for even values of and T = tt for odd A^ 
(in which case x^ = 1/A^ at r = 7r/2). 

Kitagawa and Ueda [l^ have pointed out that the non- 
linear evolution e~*^^^3 actually creates spin-squeezing, 
for r < 1/a/]V, along a particular direction. The 
maximum squeezing is obtained for the state |'0(t)) = 
e^^'^-i|V^(r)), where S{N,r) = |arctanf , A = 1 - 
(cos2r)^~^ and B = 4 sinr(cos r)^~^. We have 

[4 + (A- - 1){A - + B^)] /4(cosr)'^-^ (9) 

Equation ([9|), as a function of r, is shown in Fig. ([Ha) 
[20|. We have < i for < r < 1.15/V^ and the 
minimum, ^^-^ = l/A"^/^, is reached at r = 1.2/A'^/^. 
For <r > 7r/2, > 1 and it converges to Eq.Q, 

which eventually diverges at r = 7r/2. 

Heisenberg Limit. So far we have demonstrated that 
the non-linear evolution of a coherent spin state creates 
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FIG. 1: (color online) a) Plot of Eqs.{|7| (dot-dashed black 
line), dS} (solid blue line) and (|9]) (dashed red line) as a func- 
tion of tVN (here N = 10^). The states having x^, zeta'^ < 1 
(i.e. below the horizontal dotted line in the figure) are useful 
for quantum interferometry. b) Phase sensitivity as a function 
of the total number of particles A^t = Np. Circles are results 
of numerical simulations, the black line is the Heisenberg limit 
= 8.9/A^T, obtained for p = Popt, and the dotted blue line 
is the shot-noise = 1 / V Nt . Inset: A^ as a function of 
the number of measurements p, for fixed values of A^t- The 
optimal working point (minimum of each curve) is popt = 20, 
independently from A^t- 
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particle entanglement useful for sub shot-noise sensitiv- 
ity. This protocol has advantages when compared to the 
spin-squeezing approach discussed in [l^ for improving 
the phase sensitivity of a MZ [21]. While spin-squeezing 
is created only for short times, r < 1/a/]V, and along 
a direction 6{N, r) which strongly depends on r and 
our scheme does not require any additional rotation of 
the initial state, is fairly independent on the evolution 
time and reaches the Heisenberg limit [22[, A^hl = 1/^, 
for times for r > 1 / ^/N . Here we apply these results to 
a realistic BEC experimental setup. The coherent spin 
state can be created by splitting an initial condensate in 
two modes with the ramping of a potential barrier or by 
quickly transferring half of the particles from an initial 
condensate to two different hyperfine levels with a 7r/2 
Bragg pulse. The non-linear evolution, e~^^^ , where 
T = Ect^ Ec is the charging energy and t is the evolution 
time, is naturally provided by particle-particle interac- 
tion The non-linear dynamics of an initial separa- 
ble state has been also recently experimentally demon- 
strated with trapped ions Here we consider a Mach- 
Zehnder interferometer with input state |'^(t)) and infer 
the true value of the phase shift from the measure- 
ment of the relative number of particles at the output 
ports. These are characterized by conditional proba- 
bilities P(/i|j,^,r) = \,{j,fi\e-'^^y\i^{T))\^, being /i the 
result of a measurement. To achieve A^qcr, Eg. ([5]), 
we consider a Bayesian estimation scheme [2S, l25|.In 
Fig.([TJb) we plot the results of numerical simulations for 
r = and = 7r/2. We show A^ as a function of 

the total number of particles used in the estimation pro- 
cess Nt = Np^ being p = popt = 20 the optimal number 
of independent measurements. The circles are numeri- 
cal results (minima in the inset of ([Hb)) and the line 
is A^ = 8.9/A^T- We emphasize that the more popular 
phase estimation scheme based on the measurement of 



FIG. 2: (color online) a-c) Distributions P(/i|i, 0, r) plotted as 
a function of along circles of radius /i (taken as a continuum 
variable) at three different times during the non-linear evolu- 
tion: a) r = 0, b) T = 7r/4 and c) r = 7r/2. The solid lines 
delimit the typical size of the substructures, d-f) P(/j\j,0,T) 
as a function of 0, and for: d) r = 0, /x = 7.5 e) r = 7r/4, 
/i 2.5 and f ) r = 7r/2, /x 3.5. Here N 15. 



average moments of Jz l9[ and the corresponding error 
propagation analysis only provide shot-noise. 

Can we understand the origin of sub shot-noise without 
spin-squeezing ? Let us investigate the phase structures 
characterizing the conditional probability distributions 
P(/i|j, 6>, r), defined for discrete values of —j<fi< j. 
These distributions contain all of the available infor- 
mation about the parameter that can be extracted 
from the measurement of /i. In Figs.([2la)-([2lc) we plot 
P{jii\j,0,r), as a function of 0^ along circles of radius /i, 
at three different times during the non-linear evolution: 
©a) r = 0, ((Sb) r = 7r/4 and l^c) r = 7r/2. The 
typical size of the substructures is ~ l/>/]V in (Oa) and 
~ 1/A^ in (Ob) and (Oc) as indicated, in the figure, by 
solid lines. This is also shown in Figs.([2ld)-([2lf) where 
we plot P(/i|j, ^, r) for different /i and the same r as in 
Figs.([2la)-([2lc). The size of the relevant substructures 
indicates the smallest rotation angle needed to make the 
rotated state orthogonal to the initial one. 

Conclusion. We have explored the interplay between 
multiparticle entanglement and quantum interferometry. 
A key role is played by the quantum Fisher information. 
We obtained a sufficient condition for A/'-particles entan- 
glement, X < Eq.(|3]), which is more general than and 
incorporates the spin-squeezing condition Eq.([2]). Large 
entanglement can be obtained through a non-linear evo- 
lution and used to reach a phase sensitivity at the Heisen- 
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berg limit. Our results can have practical impact in pre- 
cision spectroscopy, atomic clock and atomic/optical in- 
terferometry and can be implemented with Bose-Einstein 
condensates and trapped ions within the current technol- 
ogy. 
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